Mean-field description of pairing effects, BKT physics, and superfluidity in 2D Bose gases 
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We derive a mean-field description for two-dimensional (2D) interacting Bose gases at arbitrary temperatures. 
We find that genuine Bose-Einstein condensation with long-range coherence only survives at zero temperature. 
At finite temperatures, many-body pairing effects included in our mean-field theory introduce a finite amplitude 
for the pairing density, which results in a finite superfluid density. We incorporate Berenzinskii-Kosterlitz- 
Thouless (BKT) physics into our model by considering the phase fluctuations of our pairing field. This then 
leads to the result that the superfluid phase is only stable below the BKT temperature due to these phase fluc- 
tuations. In the weakly interacting regime at low temperature we compare our theory to previous results from 
perturbative calculations, renormalization group calculations as well as Monte Carlo simulations. We present a 
finite-temperature phase diagram of 2D Bose gases. One signature of the finite amplitude of the pairing density 
field is a two-peak structure in the single-particle spectral function, resembling that of the pseudogap phase in 
2D attractive Fermi gases. 
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I. INTRODUCTION 

Recent experiments on two-dimensional (2D) ultra-cold 
atoms have explored many interesting phenomena includ- 
ing the Berenzinskii-Kosterlitz-Thouless (BKT) physics [1], 
superfluidity [2|, scale invariance (3), radio-frequency (RF) 
spectroscopy |4|, thermodynamics [5], pseudogap physics 
above the BKT transition temperature (6), and others. These 
experiments provide opportunities for studying more compli- 
cated 2D or layered systems related to high-temperature su- 
perconductors [7 8 1 and interface superconductivity [9|. The- 
oretical studies on dilute 2D Bose gases have been reviewed in 
Ref. [ 10 1 and those available theories are confined to weakly 
interacting regimes or temperatures close to zero or near the 
critical regime. As a consequence, finite-temperature phase 
diagrams shown in Ref. |10| are schematic instead of coming 
from a consistent theoretical description. A mean-field theory 
that works in the regime of intermediate interaction strength 
at arbitrary temperature thus would be highly desired for a 
systematic analysis of 2D interacting Bose gases. To better 
understand the physics, one needs a coherent description of 
superfluidity, the BKT transition, pairing effects, and single- 
particle excitation energy. The goal of this paper is to present 
a plausible mean-field theory with experimental consequences 
for a 2D interacting single-species Bose gas. 

For an attractive 2D two-component Fermi gas, there have 
been theories based on the phase fluctuations of the BCS the- 
ory and its extension to Bose-Einstein condensation (BEC) 
of dimers (7] QT). When the temperature T is below a pair- 
ing onset temperature, pairs with disordered phases emerge. 
When T falls below the BKT transition temperature Tbkt, 
a superfluid phase becomes stable but a genuine long-range 
ordered phase only survives at T = 0. For bosons we may 
explore similar physics. Several questions follow: How can 
BKT physics be incorporated into a theory of 2D interacting 
bosons? Does any interesting phase exist above TbktI Can 
2D bosons have an energy gap in the single-particle excita- 
tion? These issues will be addressed in a consistent theoretical 



framework. 

Ref. [10] poses a series of questions that needs to be ex- 
plored in theoretical work on 2D interacting Bose gases. The 
last one is "Could one justify a large-A*" approach which im- 
proves on existing methods by incorporating the i-matrix ap- 
proximation?". Inspired by this question, here we base our 
theory on the leading-order-auxiliary-field (LOAF) theory of 
interacting bosons [12 13 1, which is a generalization of the 
conventional large- N expansion lfl4l - [T6l . It also reduces to 
the large-N expansion in the normal phase. The LOAF theory 
is a mean-field theory for interacting Bose gases. Its advan- 
tages in describing a 3D Bose gas beyond perturbative regimes 
at arbitrary temperature are summarized below. In the follow- 
ing we will construct a mean-field theory that applies to 2D 
Bose gases and explore its thermodynamics and possible ex- 
perimental implications. 

For 3D interacting bosons the LOAF theory meets three 
important criteria by treating the pairing (anomalous) density 
field and the (normal) density field on equal footing: (i) a gap- 
less dispersion in the BEC phase, (ii) a conserving theory, and 
(iii) predicts a second-order BEC transition. Widely used the- 
ories such as the Hartree-Fock theory or the Popov theory fail 
at least one criterion 1131 . Moreover, the LOAF theory ex- 
hibits a shift in the critical temperature T c consistent with the 
results of Ref. ifTTll . We emphasize that the LOAF theory nat- 
urally recovers the Bogoliubov theory of weakly interacting 
bosons [12] and note that two Green's functions correspond- 
ing to our two density fields are indeed present in the Bo- 
goliubov theory 1 18 1. An important feature of the LOAF ap- 
proximation is that the superfluid density is closely related to 
the pairing density lfl9l and the two quantities obey a Joseph- 
son relation [20|. This will be crucial in integrating the BKT 
physics into the LOAF theory of 2D Bose gases. 

Before presenting our theory, we point out several chal- 
lenges in developing a consistent mean-field theory for 2D 
interacting Bose gases. For 2D systems with short -ranged in- 
teractions, Mermin- Wagner theorem [21] rules out the pos- 
sibility of long-range order. Therefore at finite T a genuine 



condensate is not possible in the thermodynamic limit. The 
lack of a genuine condensate implies that a gapless Gold- 
stone mode may not exist so that the dispersion of excitations 
should be gapped rather than having the gapless Bogoliubov 
dispersion relation found in three dimensions in the conden- 
sate phase. This challenge may be circumvented by confining 
the discussion on finite-size patches of the system or introduc- 
ing phase fluctuations to disorder the condensate [10]. Here 
we consider 2D Bose gases in the thermodynamic limit so we 
follow the second approach. The two-body scattering problem 
in 2D is also very different from that in 3D (seeRefs. [10, 22]). 
The scattering amplitude in 3D is finite as both the energy- 
and momentum- transfers approach zero. This allows one to 
develop an effective theory with a coupling constant set equal 
to the scattering amplitude at zero energy- and momentum- 
transfer. In contrast, the 2D scattering amplitude vanishes in 
the limit of zero energy- and momentum- transfer. Following 
the standard renormalization methods [ 15], one has to formu- 
late a running coupling constant which should be defined at a 
finite energy (or momentum) scale. A renormalization scheme 
for the running coupling constant will be presented and the 
running coupling constant indeed exhibits behavior consistent 
with the i-matrix analysis of the scattering amplitude. 

This paper is organized as follows. Section |TT| shows the 
derivation of the 2D LOAF theory. The renormalization of 
the finite temperature effective potential and the incorpora- 
tion of phase fluctuations and the BKT physics are addressed 
in detail. In Section III we present the calculated phase dia- 



gram and the experimental implications of the 2D LOAF the- 
ory. Comparisons with perturbative calculations are also pre- 
sented. Section|lV]states the conclusions of our work. 



II. LOAF THEORY OF 2D INTERACTING BOSE GASES 

The action of a homogeneous 2D Bose gas is given by S — 
J dxC, where dx = dtd 2 x and the Lagrangian density is 



C = l -W(x)H(x) + </>(x)h*<l>*(x)] - ^l^)| 4 - 



(1) 



Here h = ihdt + h 2 V 2 /2m + [i and fj, is the chemical po- 
tential. We set h = 1. A is the bare 2D repulsive coupling 
constant. After discussing its renormalization, we will show 
its connection to the 2D s-wave scattering length. This action 
is equivalent to the Hamiltonian 



2m 2 



(2) 



Although the Hamiltonian approach to BECs is more com- 
mon in the literature, having a path-integral formulation of the 
problem allows for certain expansion methods that would be 
hard to implement in the canonical formalism. The expansion 
used here is one such expansion which allows one to interpret 
a mean-field theory as the first term in a complete resumma- 
tion of the original theory by utilizing a Hubbard-Stratonovich 
transformation [23 1 and then reversing the order of integra- 
tions. Introducing auxiliary fields allows one to do the path 



integration over the original fields exactly while holding the 
auxiliary fields constant. Then one does the remaining path 
integrations by the method of steepest descent and obtains a 
loop expansion in terms of the auxiliary-field propagators. 

In utilizing the Hubbard-Stratonovich transformation we 
impose the following requirements: (i) to treat the normal and 
anamolous densities on an equal footing and (ii) to reduce in 
the weak-interaction limit to Bogoliubov's theory of interact- 
ing bosons (see [24, 25] for reviews). To insure this latter fea- 
ture we introduce the normal and pairing density composite 
fields xq an d A representing \/2Xcf>*(x)(f>(x) and \<f)(x)<f>(x) 
with the corresponding fluctuations, the Lagrangian density in 
the LOAF theory becomes lfT2l[l3l 



C = £ + [A(x)[<j>*(x)] 2 +A*(x)[<j>(x)Y 



v^Xo^Oc)! 2 



2A 



[xl{x)-\A{xt 



(3) 



Here Co denotes the kinetic energy part of Eq. ([TJ. Note that 
the pairing density field emerges from the pairing channel in 
Eq. ([TJ. These auxiliary fields introduce composite-field prop- 
agators and the counting of their loops facilitates a resumma- 
tion scheme similar to the large- N expansion |[T3~14T6ll . 

The leading-order auxiliary-field (LOAF) theory of 3D 
Bose gases has been discussed in Ref. [12, 13 1 and here we 
briefly review its derivation for a 2D Bose gas. We begin with 
the Lagrangian density after introducing the composite fields 
Xo and A as shown in Eq. (BJ. To simplify our expressions, we 
define $ = (0, <jf , \o, A, A*) T and its corresponding source 
term J = (j, j* , s, S, S*) T . The generating functional W[J] 
can be obtained from the partition function by 



sm = j 



Z[J] = e ««W> = M / v<Pe^ h ^ s W+f Jt *l 



dtd xC. 



(4) 



We note here that if we first perform the Gaussian path in- 
tegration over the auxiliary fields xo,A,A*, we recover the 
original generating functional based on the Lagrangian den- 
sity ([TJ. Instead we will reverse the orders of integrations and 
perform the (now) Gaussian integration over the fields <f>, <j>* . 
This is the gist of the Hubbard-Stratanovich transformation. 
The generating functional of one-particle irreducible (1-PI) 
graphs r[$] is the Legendre transformation of VK[J]: 



r[$] 



d d xJ\x)$(x)-W[J] 



(5) 



The action 5[$] in Eq.(p| written in terms of the auxil- 
iary fields A and xo is quadratic in the fields <j>, <fi*, so that 
the Gaussian integral in <j>, <j>* can be done exactly, leaving 
an effective action S e ff[xo,A,A*] which now depends on 
J, xo, A, A*. In order to perform the remaining integrals we 
introduce a small parameter e into the problem by replacing 
S e ff — > S e ff/e, This allows us to perform the remaining in- 
tegrals via Laplace's method or the stationary-phase approxi- 
mation. As shown in Refs. lfT3ll26l . the parameter e counts the 
number of composite-field (x and A) propagator loops. The 
LOAF approximation is the leading order in e and consists of 



just keeping the contribution at the stationary phase point. The 
details of this expansion are given in Ref . |[T3l . 

The Legendre transform of the stationary phase approxima- 
tion for Z[J] is then found to be 



T[$] = i / ' d d xd d x'ct>* a G-^{x,x')<i> h {x') 
i*Jxl-\A? '■ 



2A 



--TrllnlG- 1 ^^)]} 



Here G 1 (x, x') = S(x,x')G 1 and G l is given by 



-ihd t -h 2 W 2 /(2m) + x -A 

-A* ihd t ~h 2 V 2 /(2m)+x 



Here \ — v^Xo — M- The effective potential can be eval- 
uated by a Wick rotation to the imaginary time r — > it 
and using the standard Matsubara frequency summation. For 
homogeneous static fields we define the effective potential 
V e ff = r[$]/V/3, where V is the volume of the system and 
/3 = l/(ksT). The inverse Green's function becomes 



G~ 



-iu n + Cfe + X 

-A* 



-A 



Here uj n is a bosonic Matsubara frequency and efc = 
h 2 k 2 /(2m). The last term in V eff is iTr{lnG -1 } = 

i£ fc ££>K + " f = SJK/2 ) + (V/?)in(i - 

e _,3 " fc )]. Here w fe = V( e *; + x) 2 ~ l^l 2 - Then one obtains 
the expression for the (unrenormalized) effective potential 



Veff = X\<f>\ 2 



2 



Wfc 

2 



iMl 



2 



(x + /j) 2 

4A 



2A 



(6) 



At the minimum of the potential 8V e ff/S<j>* = which 
leads to the condition x4> ~ A(f>* = 0. Thus if the minimum 
is not at cp = 0, we can use the U(l) symmetry to make <f> 
and A real and then we obtain in that case \ — A. Indeed, 
at T = the minimum is not at zero and the single-particle 
Bose-Einstein condensation (BEC) corresponds to a finite (f> = 
•v/po. Therefore in the presence of the single-particle BEC, 

X = A and the dispersion is gapless u>k = \/ efc ( efc + ^X)> 
where x plays the role of \po in the Bogoliubov theory of 
weakly interacting bosons. 



A. Renormalization of LOAF theory 

So far the theory has not been renormalized and it is ultravi- 
olet divergent. In Popov's approach to 2D Bose gases 1 10, 27 1 
a high-energy cutoff has been introduced and this cutoff re- 
mains in the equations of state. Here we systematically renor- 
malize the coupling constant, chemical potential, and vacuum 
energy so that there is no ultra-violet divergence in our the- 
ory. By examining V e ff shown in Eq. (|6]l one can see that the 



ultraviolet divergence comes from the integral of Wfe. The ul- 
traviolet divergent part of w^ at large k is efc + x — \A\ 2 /(2ek) 
in both 3D and 2D. In 3 D one can relate these terms to usual 
renormalizations evaluated at fi = q = 0, where Ml and q 
are the energy transfer and momentum transfer of two-body 
scattering. In two dimensions, however, one needs to de- 
fine the renormalized coupling constant at a finite value of fi 
and/or q since the scattering amplitude vanishes at zero ft and 

qMM- 

In our resummed theory the renormalized running coupling 
constant is defined as the inverse of the second derivative of 
the effective action with respect to the field A. This leads to 
an equation 

X^KT) =2DAUq ^ n (?) 

The inverse AA* propagator can be obtained by 

6 2 r[<S>}/SA5A*. Explicitly, 

°aa^^^)- 2A +2SA(x)6A*(x'y (8) 
After a Fourier transform this becomes 

D^.( g ,n,r)=^-in(g,n,T). (9) 

Here H(q, O, T) is the bubble diagram discussed in Ref. |fl9l . 
The divergent part of U(q,fl = 0,T = 0) in 2D is 
(1/2tt) / dkk[y/e 2 + {e q /2) 2 ]- 1 . We define the physical 
renormalized coupling constant in 2D as the running coupling 
constant at T = 0, O = 0, q = qo (of course we can choose 
any other scale to define the renormalized coupling constant). 
For a comparison with the relativistic case, one can refer to 
Ref. [28 1 and references therein. 

We define an intermediate renormalized coupling constant 
by just keeping the ultraviolet divergent part of the bubble di- 
agram. 



dkk 



X R (q ) A J 2k ^/ e 2 + (e 9o /2)2 



(10) 



In terms of this intermediate renormalized coupling constant 
(which we will choose to be our definition of the renormalized 
coupling constant A# ) one finds that the running coupling 
constant at any scale is given by 

1 ' (n(g,n,T)- 



X R (q,n,T) X R (q ) 
dkk 



1 



^ V4 + (^m 2 , 



(id 



In three dimensions one is allowed to set qo = 0, since the 
two-body scattering is finite at zero momentum transfer. 

The connection between the 2D renormalized running cou- 
pling constant and the t-matrix calculations of the 2D scat- 
tering amplitude can be seen clearly at T = 0. One can 
show that U(q,n = 0,T = 0) - H(q ,fl = 0,T = 0) = 
(l/27r) \n(q/qo) and as a consequence 

Ah (g ) 



X R (q,n = 0,T = 0) 



1 



A- In (A) 



(12) 



Therefore X R (q,n = 0,T = 0) ->• as q ->• 0, which is 
consistent with the finding that the 2D scattering amplitude 
vanishes at zero energy- and momentum- transfer ifTOl l22l in 
the i-matrix calculations. This analysis justifies our choice of 
a finite momentum transfer in defining the running coupling 
constant. The 2D scattering amplitude at a small momentum 
transfer is given by f 2D ~ 47r/[rn| ln(pa 2 )|], where a is the 
2D s-wave scattering length flOl |22l . After the coupling- 
constant renormalization, A.r may be set to be proportional 
to f 2D and we formally write 



A, 



Anh 2 



m\ In (pa 2 



(13) 



Here a is the scattering length measured at the momentum 
transfer q an d T — 0. We focus on the regime pa 2 < 1, 
which is consistent with the requirement of a dilute gas. This 
less restrictive condition, when compared to the condition 
| In | ln(pa 2 ) 1 1 <C 1 for those perturbative methods reviewed in 
Ref. lfTOl . allows us to explore intermediate-coupling regimes 
where lln^a 2 )!" 1 ~0(1). 

In contrast to the renormalization of the coupling constant, 
the chemical-potential and vacuum-energy renormalizations 
are found to have the same form in two and three dimensions. 
The chemical-potential renormalization is 

PR _ /f 

Xr A 



k 



1. 



(14) 



The vacuum-energy renormalization is given by 



v (i2r -v r 

Vr -4^- V °-4\ ^2 

■ k 



E 



Here Vq and Vr are the bare and renormalized vacuum energy. 
The LOAF theory written in terms of the renormalized pa- 
rameters is thus given by 

A*£ _ A(cf>*) 2 _ ( X + /j) 2 W 

2A 



Veff = x\^ 



1 

27 

1 



dkk 



Wfc 



2 

Cfc 



4A 



x- 



\A\ 2 



V £fc+ K>/ 2 ) 2 



+ 



ln(l 



o~P"k 



(16) 

Here we drop the index R and the vacuum energy and use 
continuity to extend our renormalization procedure to the re- 
gion where \ 7^ A. The equations of state are given by 
SV ef f/Sx = 0, 5Veff/6A* = 0, and -SV e ff/Sn = p, which 
are explicitly shown in the main text. 

The LOAF theory then only keeps the leading order terms 
and one obtains the equations of state (EOS) by minimization 
of the renormalized 14//: 



.4 



= Pa + A 



d 2 k 



l + 2n(uj k ) 
2uj k 



1 



P = Po 



d 2 k 



e ^i + M« k )}-- 2 



V4 + (s/2) 2 
i" 



(17) 
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Figure 1. (a) The condensate fraction as a function of | ln(pa 2 )| _1 
at T = (solid line). The dashed line shows the prediction from 
Eq. l |19fr . The inset shows the pairing density field A as a func- 
tion of j In (pa 2 ) | _1 . (b) The chemical potential as a function of 
ln(pa 2 )| _1 at T = (solid line). The dashed line shows the pre- 



/i<->. diction from Eq. (T8D . 



Here \ = V^Xo - M. w fe = V( e k + X) 2 - \A\ 2 , n(x) = 
\j{er x — 1) is the Bose distribution function, and p = 4% 
is the condensate density, where <f>Q denotes the expectation 
value of <p. The density is related to the chemical poten- 
tial via p = (x + p)/2X. In addition to the EOS, there 
is a BEC condition \4> ~ Ac/)* — 0. We define k 2 = p, 
Eq = h 2 k 2 /2m = k R To and use ko, Eq, and To as our units 
with fog = 1. The quantity qo is the scale that one defines 
the running coupling constant X(q, il,T) in a 2D Bose gas. 
It also defines the momentum transfer at which the scatter- 
ing length a is measured. We emphasize again that in 2D 
qo cannot be chosen to be zero due to the peculiar behavior 
of the scattering amplitude. However, this choice of scale 
does not effect any physically measurable quantities as argued 
in standard renormalization methods [15J. Here we choose 
<7o/&o = 0.1. To compare with experimental parameters, one 
only needs to find out the momentum transfer q at which the 
relation f 2D ~ 47r/[m| ln(pa 2 )|] is defined. That qo then 
serves as the renormalization energy scale in the EOS. 

At T = the LOAF theory predicts a finite p Q and BEC. 
Moreover, the BEC condition \(\> — A<\>* = requires that 
X = \A\ so that w fc = v /e fe (e fc + 2\) is gapless. This gap- 
less excitation is associated with the Goldstone mode in the 
BEC phase [19). The equations of state become (A/X) = 
Po - (^/87r)ln(2A/e go ) and p = p + A/8n. Figure [j] 
(a) shows that po/p decreases and A/ksTo increases as the 



interaction | ln(pa 2 )| _1 increases. The chemical potential is 
shown in Figure[TJ(b). The condensate fraction and the chem- 
ical potential from perturbative calculations [29 1 are given by 



/' 
Po 



AiTh 2 p 
m\ ln(pa 2 )| 
1 



[l + 0(\Hpa 2 )\-% 



P |ln(pa 2 )l 



0(\Hpa 2 )\- 



(18) 



(19) 



The dashed lines in Figure [TJ show the results from pertur- 
bative calculations. One can see that for extremely small 
| ln(pa 2 )| _1 , the LOAF theory agrees well with the pertur- 
bative calculations. Thus the LOAF theory agrees quantita- 
tively with perturbation theory results and should be reliable 
in the weakly interacting regime. Importantly, the LOAF the- 
ory could help explore interesting physics beyond the weakly 
interacting regime. 



B. Phase fluctuations 

At finite T, Mermin-Wagner theorem |21 1 rules out the pos- 
sibility of long-range orders so BEC cannot survive. In other 
words, the U(l) symmetry of the Lagrangian densityfTlshould 
remain unbroken at any finite T. The lack of a broken sym- 
metry implies that a gapless Goldstone mode should not exist 
and the excitation should be gapped. This is consistent with 
Eq. ( fT7| ) since the equations cannot be satisfied by a gapless 
dispersion. In the LOAF theory this indicates that the BEC 
condition cannot be met so po must vanish. Eq. ( fT7j , how- 
ever, does not rule out the possibility of a finite A and indeed 
we found finite values of A at finite T, The finite expectation 
value of A implies a diatomic condensate, which also could 
break the U(l) symmetry of the Lagrangian density ([TJ [19| 
and violate the Mermin-Wagner theorem. 

This apparent conflict was analyzed in great detail by Wit- 
ten 113011 in the SU(N) Thirring model at large- N. His argu- 
ments would lead here to the conclusion that there is a gapless 
mode associated with the phase of A but that this will not be a 
Goldstone mode. He also argues that apart from the fact that 
there is no symmetry breaking at finite T, the fermions (which 
in our case will be the boson </>) will have an energy gap in the 
phase-disordered regime and that apart from the details of the 
BKT phase transition [31], the large-A*" predictions (which are 
similar in spirit to the LOAF theory) are quite reliable. 

In the BCS theory of attractive fermions, the method of ob- 
taining a BKT transition in 2D by introducing phase fluctua- 
tions to disorder the would-be order parameter has been dis- 
cussed in Refs. |7J[TTJ. There after obtaining the BCS the- 
ory at the mean-field level, one further introduces a phase for 
the gap which is the analogue of the anomalous condensate 
present in the LOAF theory. It has been shown that the BCS 
theory can be derived from the same LOAF framework 11321 
and one can obtain similar Josephson relations in both the 
BCS theory of attractive fermions and the LOAF theory of 
repulsive bosons ll20l . Therefore phase fluctuations may be 
incorporated into the LOAF theory of bosons using a method 
similar to that for the BCS theory of attractive fermions. 



Following the phase-fluctuation method of the Thirring 
model and the BCS theory, we introduce phase fluctuations 
into the solution of the EOS of the LOAF theory 1171 FITl |30l . 
This procedure also introduces the BKT transition to our the- 
ory and determines where the superfluid phase is stable. The 
idea is to include a fluctuating phase in the pairing density 
field so it becomes Ae 10 . As discussed in Ref. |20|, the phase 
of the pairing field is twice of the phase of the bosonic field 
and both phases originate from the U(l) symmetry of the 
Lagrangian density ([TJ. The amplitude A is determined by 
the EOS, Eq. {FT), and can be finite. Similar to Ref. Q 
we work with the minimum where A is constant. Follow- 
ing Refs. QUO), in the action containing 9 we only keep the 
leading-order contribution of the phase fluctuation, which is 
proportional to J d 2 x(V6) 2 . This kinetic -energy term of 6 
corresponds to the XY model and thus the conventional BKT 
physics On is introduced to the theory following the argu- 
ments shown in Section 5 of Ref. |7|. The proportionality co- 
efficient of the kinetic -energy term of 9 is the phase stiffness, 
which is equal to the superfluid density p s = 5 2 V e ff/5vSv 
following a similar calculation of Ref. [20|. Here v — V0. 

Similar to the results presented in Ref. 0, the phase fluctu- 
ations obey (e 10 ) = and as a consequence, there is no long- 
range coherence of the pairing density field. The correlation 
of the phase is (e i6( - x h i6 ^) _ | x |- mT /( 2 ^) below the BKT 



i9(x) p i8(0)\ 



-|x|/x 



transition temperature Tbkt and (e 
above Tbkt, where xo is the characteristic length for the de- 
cay of the correlation (ST). The conventional picture is that 
vortex-antivortex pairs are bound below Tbkt and a super- 
fluid phase is possible. Above Tbkt vortices and antivortices 
unbind and proliferate so the superfluid is destroyed. Thus 
the phase fluctuations introduce a BKT transition separating a 
low-T superfluid phase and a non-superfluid phase. 



III. RESULTS AND DISCUSSIONS 
A. Finite-T phase diagram 

We now construct the finite-T phase diagram. The pairing 
onset temperature T p is determined by the EOS ( fTTJ when the 
amplitude A first becomes finite. Figure [2] (a) shows T p as a 
function of r\. There is no genuine phase transition across T p . 
Below T p bosons form composite pairing density field but the 
phase is random. In other words, below T p there is a phase- 
disordered diatomic quasi-condensate. We note that in the 
2 + 1-dimensional relativistic O(N) <p A model, there is a reso- 
nance in the <fi — (j) scattering amplitude near twice the mass of 
the scalar meson at small coupling [28]. Further studies may 
determine if there is finite two-body binding energy. Fig. [21(b) 
shows the growth of the amplitude of A as T decreases. 

Next we investigate where superfluidity becomes stable. 
The method for determining the superfluid transition temper- 
ature comes from the analogy with the XY model once one 
adds the phase fluctuations 0. According to the theory of 
BKT transition [31 1, the superfluid is unstable above Tbkt 
due to vortex-antivortex proliferation, where Tbkt is deter- 
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Figure 2. (a) Phase diagram of a 2D Bose gas from the LOAF theory 
with phase fluctuations. The solid and dashed lines show the BKT 
and pairing onset temperatures. Regions I, II, and III correspond to 
the normal, pairing, and superfluid phases. The details of physical 
quantities along the dotted line (j ln(pa 2 )| _1 = 0.5) are shown in 
panel (b). Inset: The gap in the dispersion, u)k(k = 0) = uiq, as 
a function of T for | ln(pa 2 )j _1 = 0.5. (b) The amplitude of the 
pairing density field A and the superfluid fraction p s j p as a function 
of T for | ln(pa 2 )| _1 = 0.5. The dashed line is (2m/7r)T/T and 
when it intersects the curve of p s /p, the BKT transition occurs and 
the superfluid density drops to zero, (c) Comparison of (1) T^ from 
Eq. ^ , (2) T p from the LOAF theory, (3) Tbkt from Eq. ((23) with 
C = 2.152, and (4) Tbkt from the LOAF theory in the weakly 
interacting regime. 



mined by 



^bTbkt — 7. — Ps{Tbkt)- 
Zm 



(20) 



The superfluid density p s of the LOAF theory has been dis- 
cussed in Ref. [19] and it can be obtained from the Landau 
two-fluid model or the current-current response function. It 
has been argued that the pairing density field is crucial in sus- 
taining a finite p s . This feature is similar to the fermionic 
BCS theory, where the superfluidity comes from Cooper pairs. 



Ps= P 



/i 2 

rn 



d 2 k 



dn(u> k ) 
duj k 



(21) 



Fig. 0(b) shows Ps/ p as a function of T for | ln(pa 2 )| = 0.5. 
When this curve intersects the line of (2m/h 2 Tr)T/To, the 
BKT condition p0| ) is met and vortex- antivortex proliferation 
will destroy the superfluidity above Tbkt- As a consequence, 
p s jumps to zero and the superfluid phase is only stable below 
Tbkt- 

Following this procedure we show Tbkt as a function of 
| ln(/9a 2 )| _1 in Fig. 0(a). There is a genuine phase transi- 
tion across Tbkt since the superfluid density is discontinu- 
ous across this boundary. We therefore identify three different 
phases of a 2D Bose gas at finite T as shown on Fig. [2] (a): 
Regime I above T p corresponds to a normal gas with no pair- 
ing density nor superfluidity. Regime II in between Tbkt 
and T p is a non-superfluid phase with a finite amplitude of the 
pairing density field but no phase coherence. Regime III be- 
low Tbkt is a superfluid phase with algebraically decaying 
phase correlations. One can see that T p is about 2 — 3 times of 
Tbkt and is not a high-temperature scale as one may inferred 
from Ref. [ 1 ]. In the LOAF theory, the phase diagrams at fixed 
values of | ln(pa 2 )| _1 are similar to the one shown in Fig. 
(b) (i. e., the amplitude of A smoothly approaches zero at T p 
while p s drops to zero at Tbkt) but with different values of 
T p and Tbkt which can be inferred from Fig. [2] (a). 

The superfluid transition temperature in the weakly inter- 
acting regime has been theoretically studied in previous work 
(see ifTUl for a review). In Ref. |[33l a renormalization-group 
approach predicts that 



T 



RG 



2irp 
mln | ln(pa 2 



(22) 



Monte Carlo simulations in Refs. J34l l35ll suggest the func- 
tional form 



T, 



BKT 



2irp 



1 



m C + ln|ln(pa 2 )|' 



(23) 



Fitting Tbkt from the LOAF theory at weak coupling to the 
functional form of Eq. ( |23| l leads to the results shown in Fig- 
ure [2] From the fitting we obtain C = 2.152 as compared 
with the Monte Carlo result of Ref. ||35l , which yielded the 
numerical value C = 3.409 ll36ll . Figure [27c) shows T p and 
Tbkt from the LOAF theory along with T^ G and T B kt 
from Eq. ( f2"3"| > with a constant C = 2.152 in the weakly in- 
teracting regime. One can see that the LOAF theory predicts 
a Tbkt that agrees well with the functional form suggested 
by Monte Carlo simulations when |ln(/3a 2 )| _1 < 0.01. In 
contrast, T^ G seems to agree with T p when the interaction 
strength is extremely small. Moreover, T p and Tbkt from the 
LOAF theory approach each other as | ln(pa 2 )| _1 approaches 
zero. As the interaction becomes stronger, Tbkt as obtained 
from the LOAF theory starts to deviate from Tbkt from 
Eq. (|23j and remains smaller than both T c and Eq. ( |23] l. 
The double logarithmic behavior of both T p and Tbkt from 



the LOAF theory in the weakly interacting regime also sug- 
gests that our theory captures the qualitative features of both 
Monte Carlo and renormalization-group results. LOAF theory 
serves as an extrapolation beyond weak coupling and provides 
predictions that can be compared with experiments. 

Interestingly, the phase diagram of Fig.[2](a) is similar to the 
phase diagram of a 2D Fermi gas with attractive interactions 
ifTTI . There are subtle differences [33. For example, Tbkt 
for fermions increases as the attractive interactions increase 
but for bosons it increases as the repulsion increases. The slow 
increase of Tqkt as a function of |ln(|0a 2 )| _1 away from 
the weakly interacting regime agrees qualitatively with more 
recent Monte Carlo simulations [38 1. 



B. Other experimental implications 

We now address the issue whether a finite amplitude of the 
pairing density field results in any observable effects. One 
observable signature is a gapped excitation energy spectrum 
Uk = \/(e& + x) 2 ~ l^l 2 at finite T (consistent with the anal- 
ysis of Ref. |[30"1 ), which is different from the gapless Bo- 
goliubov spectrum. The inset of Fig. [2] (a) shows the gap 
uJk(k = 0), which is finite and increasing as T increases. The 
dispersion may be measured using Bragg scattering [39|. An- 
other possible signature may be revealed by the analogue of 
radio-frequency (RF) spectroscopy, which shows the potential 
of measuring the spectral function, or equivalently, the imag- 
inary part of the single-particle Green's function [40 1. For 
fermions with attractive interactions, the spectral function of 
a homogeneous gas could show a two-peak structure due to 
the particle-hole mixing in the formation of Cooper pairs [40] . 
Here we investigate if there is a similar structure for bosons. 

The single-particle Green's function from the LOAF theory 
is Q9) 



G MF (k,ioj n ) 



iu n + e k + x 



(24) 



Here u n is the bosonic Matsubara frequency. Making the an- 
alytic continuation iu n — > uj + i0 + , one obtains GMF(k,ui) 
ifTHl . The spectral function is defined as A(k,oj) — 
2haG(k,ui), where G(k, ui) is the single-particle Green's 
function of a given theory, and it satisfies the sum rule 



dtO 

2^ 



A(k,u) = 1 



(25) 



for any T. 

The LOAF theory predicts a qualitative difference in the 
spectral function when the amplitude of A becomes finite. 
Below T p , one has Lu k = \J{e k + X + ^)( £ fc + X ~ 3) an d 
A(k, u) = 2iru 2 k 5(uj — Wfe) — 2ttvIS(uj + uj k ). Here u%\, v\ = 
[( e fe + x)/ UJ k ± l]/2- This expression implies that there are 
two peaks at ui = ±u k for a fixed k. This is in contrast to the 
spectral function in regime I where A = 0. In regime I, oj k — 
€ k + X so the spectral function is A(k,u>) = 2ttS(oj — Wfc). 
There is only one peak in the spectral function in the normal 
phase. 



The different numbers of peaks in the spectral function be- 
low and above T p is a prediction from the LOAF theory. The 
delta-function peaks are due to the fact that only the leading- 
order contributions from the composite fields are included at 
the mean-field level. To obtain the widths of the spectral 
peaks, one has to go beyond the leading-order theory. For the 
relativistic <f> A theory it has been shown in Ref. [41 1 that by 
considering a self-consistent e expansion for the generating 
functional of the two-particle irreducible graphs, one obtains 
an approximation to the coupled Green fuinction equations 
which generate a finite imaginary part in the self energy and 
as a consequence the spectral peaks will be broadened. 

In particular, above T p when A = 0, our theory beyond the 
leading order can be formulated as a set of Schwinger-Dyson 
(SD) equations following Ref. BTll 

G~ 



G MF 



GVA. 



V- 1 = V~ n 



n. 



n = / GGK. 



(26) 



Here S is the self energy of bosons, V is the composite- 
field propagator, A is the interaction vertex, Dq 1 (x,x / ) — 
6(x — x')/\ is the inverse bare composite-field propagator, 
n is the self energy of composite fields. To lowest order one 
may use the approximation A = 1, T> = T>q, G = Gmf, a nd 
II = no = J GmfGmf on the right-hand side. By writing 
£ = ReS + ilmE, one can see that the delta function in the 
spectral function is replaced by a Lorentzian function with a 
full width at half maximum equal to 2Im£. Below T p there are 
two finite composite fields, A and x, a nd the Green's functions 
of (f>, x, a nd A are mixed to form a 5 x 5 matrix G a p, where 
a, j3 = 1, • • • ,5 corresponding to (cf>, </>*,%, A, A*). Follow- 
ing Ref. ATI , the SD equations are 



G 



a/3 



^MF,a/3 + ^ap, 



^a/3 — J ^aa'P 



•G a i a iiGwRnh. a 



x"/3"/3, 



A, 



a/37 



*a/37 



S^al: 



(27) 



Here repeated indices are summed, A Q/ g 7 is the bare vertex 



function, and 4> 1 is an element of ( 



,X,A,A*). One may 



implement further approximations such as the bare-vertex ap- 
proximation and resort to numerical methods for evaluating 
the Green's function beyond the mean-field level as discussed 
in Ref. EI]. 

Instead of performing these intricate numerical calculations 
that should only lead to quantitative corrections, here we fo- 
cus on the qualitative features of the spectral function already 
present in the mean-field level but modified by the broadening 
of the higher-order corrections. The delta functions in A(k, uS) 
evaluated from Gmf will be broadened by effects beyond the 
LOAF theory and for illustrative purposes only, we introduce 
a Lorentzian function by 8{x) — > („) x 2 + / r / 2 )2 to approxi- 
mate the spectral peaks. The width V could be obtained from 
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Figure 3. Spectral functions at fixed k - k for (a) T > T p (T/T = 
8), (b) T P >T> T B kt (T/To = 4), and (c) T < T BK t (T/T = 
1.3). They belong to regimes I, II, and III of Fig. [21(a), respectively. 
Here|ln(pa 2 )|" 1 =0.5. 



the full expression of the self energy £ outlined above, but 
here for simplicity we set T/E = 0.02. A fully numerical 
calculations similar to the T — calculations in Ref. Il42l 
may help determine the dependence of T on T and the inter- 
action strength. Our approximation of the delta function still 
respects the sum rule ( f25| > for A(k, w). Figure [3] shows the 
spectral function A(k,u>) at fixed k = ko for three selected 
temperatures corresponding to regimes I, II, and III of Fig. [2] 

When T > T p , A — so there is only one peak as shown in 
Fig. 151 a). Below T p the finite amplitude of the pairing density 
field induces another peak in the negative energy region. In the 
Bogoliubov theory of weakly interacting bosons, the Bogoli- 
ubov transformation mixes the creation and annihilation op- 
erators. Since the LOAF theory is a natural generalization of 
the Bogoliubov theory, the pairing density field includes sim- 
ilar mixing effects. Thus the spectral weights of positive and 
negative energy states are correlated. The two-peak structure 
reflects this type of correlation effects. The spectral function 
of a T — 2D Bose gas has been evaluated using numeri- 
cal functional renormalization group method focusing on the 
positive-energy peak l42l . Since the single-particle Green's 
function |24) only contains information about the amplitude 
of the pairing density field, it does not exhibit observable sig- 
natures of the BKT transition (see Fig.[3j(b) and (c)). 

We note that the spectral function of a fermionic BCS su- 
perfluid is A(k, u) — ul5(uj - E k ) + v\5{w + E k ), where 



[1 ± (c fc - li)/E k ]/2, E k = y/(e k - m) 2 + A 2 , and 
A is the gap function. When A > 0, there are two positive 
peaks reflecting the pairing between fermions. The spin statis- 
tics, nevertheless, causes one positive and one negative peaks 
for bosons. We emphasize that although this negative peak 
should also survive in 3D Bose gases l43l . its appearance in 
2D Bose gases is a more direct evidence of the many-body 
pairing effect because the BEC vanishes at finite T. 

Spectroscopies probing single-particle excitations such as 
RF measurements may be sensitive only to the existence of 
an energy gap but not to phase coherence [40 1. From Fig. [2] 
(b) and Fig. |3jwe reach a similar conclusion. To probe the 
BKT transition and the superfluid phase below it, in addition 
to Ref. [ 1 1 we suggest experiments such as the measurement 
of the second sound, which has been shown to be an indication 
of superfluidity using hydrodynamic approaches in both 3D 
fi4ll and quasi- ID |45l geometries and should have the same 
resolution in 2D. 

IV. CONCLUSION 

In summary, we present a coherent mean-field picture of 
pairing effects, superfluidity, BKT physics, and single-particle 
excitations by introducing phase fluctuations into the LOAF 
theory of a 2D interacting Bose gas in a manner similar to 
what is done in the 2D BCS theory as well as in the Thirring 
model. In addition to mapping out the phase diagram at finite 
T, our theory predicts observable signatures of pairing effects 
above the BKT transition temperature, which resembles the 
pseudogap physics of 2D Fermi gases (6). The LOAF the- 
ory agrees reasonably with previous results from perturbative 
calculations, renormalization-group calculations, and Monte 
Carlo simulations in the weakly interacting regime but fur- 
ther explores the regime of intermediate interaction strength. 
By implementing the local density approximation for trapped 
gases, our theory may provide more insights into experiments 
suchasRefs. (3JI3- 
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